Three-dimensional elastic reverse-time migration using a high-order temporal and spatial staggered-grid finite-difference scheme

Three-dimensional (3D) elastic reverse-time migration (ERTM) can image the subsurface 3D seismic structures, and it is an important tool for Earth’s interior imaging. A common simulation kernel used in ERTM current staggered-grid finite-difference (SGFD) method of first-order wave equation. However, mere second-order accuracy time SGFD bring non-negligible dispersion, which reduces further leads to distortion imaging results. This paper proposes a vector-based using high-order obtain high-accuracy images. approach new high-resolution workflow that improves conventional from numerical simulation. The proposed established on quasi-stress–velocity equation its vector wavefield decomposition form. Advanced schemes their corresponding coefficients with fourth-order temporal solve quasi-linear system. normalized dot product condition produces high-quality images wavefields solved method. Through examples, we test efficiency analyze how simulations affects quality. We include obtains highly accurate